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Abstract 

The Fade approximant technique and the variational Monte Carlo method 
are applied to determine the ground-state energy of a finite number of charged 
bosons in two dimensions confined by a parabolic trap. The particles interact 
repulsively through a Coulombic, 1/r, potential. Analytic expressions for the 
ground-state energy are obtained. The convergence of the Fade sequence and 
comparison with the Monte Carlo results show that the error of the Fade 
estimate is less than 4 % at any boson density and is exact in the extreme 
situations of very dilute and high density. 
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The interest in finite boson systems which are confined in traps increased suddenly af- 
ter the experimental observation of signals of Bose-Einstein (BE) condensation in clusters 
of alkali-vapour atomsS These are confined three-dimensional dilute systems interacting 
through short-range potentials, which basically can be described by a modification of Bo- 
goliubov's theory.H 

Lower- dimensional systems as well as bosons interacting through long-range potentials 
are also widely studied. These studies are motivated, for example, by proposals of obtaining 
quasi two-dimensional clusters by making the traps highly anisotropic! Theoretically some 
regularities of the spectra of energy levels have been identified to check experimentally the 
effective quasi-two-dimensionality of the system.! A particularly interesting aspect of trapped 
systems is related to the fact that quantum fiuctuations do not destroy BE condensation 
in lower dimensions.! On the other hand, the ground-state properties of the homogeneous 
charged Bose gas has been studied in 3D.Q Vortex systems in high-Tc superconductors have 
been shown to be equivalent to 2D bosons with logarithmic interactions.0 

In the present paper, we study a model of two-dimensional charged bosons confined in a 
parabolic trap. Instead of an often used long-wavelength field theory description, we start 



from the A^-particle hamiltonian, which in harmonic oscillator units (i.e. y/i/(mc(Jo) for 
length, hujQ for energy, etc.) reads 

TT -I N 1 



where ujq is the dot confinement frequency, (3 = y {^^) / {huo) , and m and q are, respec- 
tively, the particle mass and charge. By varying (3 we modify the "density" of the system. 

Following our approach for electrons in a parabolic quantum dot,! we will construct two- 
point Fade approximants for the ground-state energy of the Hamiltonian (Q) from the weak 
and strong interaction expansions of the energy 

el^^o =bo + hf3'' + Kf3' + ..., (2) 
el^^oo =/5'{ao + a2//3^ + ...}. (3) 

When jS << 1, we have a system of weakly interacting bosons in a harmonic potential. 
For (3 = we have non-interacting bosons with energy bo = N which is the leading 
contribution to the energy (|^). Ferturbation theory is applied to compute 63 and bg, yielding 

,3 . < OIV^IO >^ ^^/|. (4) 

^ _ ^ < 0\V\int >2 

= ciN{N -If + C2N{N -I), (5) 
ci = -0.05300695, cs = -0.11982050, 

where V = J2i<j l^i ""O l"^' ^ shorthand notation for the harmonic oscillator inter- 

mediate states respecting conservation of angular momentum (equal to zero in the ground 
state). Due to the pair character of the potential, only one- and two-particle excitations 
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are allowed, eoiint) is the energy of the non- interacting intermediate state and eo(0) = A^. 
We would like to notice that in order to apply the Pade approximant technique one needs 
to calculate the coefficients with rather high accuracy. This was achieved in the following 
way. The one-particle excitation part (namely the coefficient Ci) was calculated from the 
straightforward summation of a single sum which converges rather fast. The convergence of 
the double sum in the two-particle excitation part, however, is slow. Therefore we used an 
alternative approach and calculated the coefficient C2 from the solution of the two particle 
interaction problem which is given in the Appendix. 

On the other hand, when (3 — > oo, a scaling of coordinates r — > /5r in the Hamiltonian (1) 
showsi^ that the potential energy (Coulomb repulsion plus parabolic confinement) behaves 
as whereas the kinetic energy ("fluctuations") is of order zero in (3. Consequently, in this 
limit the system behaves as classical point particles and oq is the minimum of the classical 
potential energy. The next term 02 takes account of the zero-point fluctuations, which is 
given by 02 = Y.a^al'^- The coefficients Oq and 02 for A^ = 2 — 5 are given in Ref. [1^. 
They were obtained from a strong-coupling expansion of the exact quantum-mechanical 
Hamiltonian. For larger A^, ao and a2 were obtained by minimising the potential energy 
with a combination of Monte Carlo and Newton methods and the normal frequencies were 
computed by solving a classical small-oscillation problem.0lli Fitted expressions (as given 
in Ref. |^) for the coefficients and 02 for 6 < A^ < 210 are 

aQ/N^I^ = 1.062 - 0.875/ArV2 _ g.lSS/AT, (6) 
a2/A^^/^ = 0.573 + 0A75/N^/^ - 0.160/Ar. (7) 

These expressions also have the correct N ^ 00 asymptotic behavior. 

Two-point Fade approximants {Ps,t{P)} are used to estimate the ground-state energy. 
They are quotients of polynomials reproducing the first s -|- 1 terms of the expansion (^ 
and the ffist t + 1 terms of the (3 ^ 00 series (^). By construction, the approximants are 
asymptotically exact in both the (3 —>■ and (3 ^ 00 limits. As it is common with the Fade 
approximants, we shall check the convergence of a sequence {Ps^t{f3)} running parallel to the 
diagonal {Ps^s(/5)}.lll In the present problem we use the off-diagonal sequence {PK+3,K{f3)}, 
which is the ffist non-trivial sequence with no singularities. We will show that it exhibits 
good convergence properties at intermediate (3 values. The explicit form of the ffist elements 
of this sequence are the following 

P3,o(/3) = + . ^y' , (8) 

1 + hf3/ao 

PaM = bo + { 1 - ^-^^f—jTTT.. \ ' (9) 
qi = bs/ao, 

P.m = bo + bsf3' \ 1 - ""^2^ \ ' (10) 

l + qil3 + q2(3^ + qsp-^ . 

qi = 0092/^3, 93 = ^392/00, 

'02 - 60 
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ao 
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(02 - bo)\bQ\/ao + bj/ao 



qs 



al/bs + 263(02 - bo)/ao' 



In Figs. |l|(a) and |l|(b) we show the relative differences between consecutive approximants 
for 20 and 210 bosons, respectively. From these figures, the maximum error of the P6,3(/3) 
approximant is estimated to be lower than 4% for any value of f3. For N = 20(210) the 
maximum error is 3(4)% which is reached for (3 = 0.8(0.4). 

The explicit form of the coefficients and bk suggests that e/N is an "almost universal" 
function of the variable N(3^ for large A^. We found that for > 90, the coefficients reach 
their asymptotic forms, and we have 

bkP' ^ Nh{NP^f'\ (12) 
au/P^ ^ N^'^hk/{Np^fl^, (13) 

where 60 = 1, &3 = i^f , h = -0.05300695, fio = 1-062 and fis = 0.573A^i/^. The coefficients 

bk and are numbers, except for the coefficient 0,2 which exhibits a smooth dependence on 
A^. In the interval 90 < < 210 02 is approximately a constant (^ 2) and consequently for 
the energy we obtain 

e{P)^Nf{Np'), (14) 

where / is only a function of Np^. The P6,3(/3) estimates for A^ = 20, A^ = 90 and A^ = 210 
are given in Fig. |^, showing the "approximate" scaling when A^ > 90. Note that bosons 
interacting through short-range potentials also show scaling behaviour in the large- A^ limit. 
Indeed, in this limit thev are described by the Gross-Pitaevskii equation, which may be 
written in a scaled form.Ej 

In order to have an independent check on the accuracy of the present Fade approximant, 
we carried out variational Monte Carlo (VMC) calculationsEl for the ground-state energy. 
The trial wave function used in the computations was (up to a normalisation constant) 




*T= n</'(^o (15) 



where 0(r) = e ''^/^ is the ground-state function of a boson in the harmonic potential, and 
the pseudopotential u{rij) was chosen as 

uir) = -(A — r r^ln r), (16) 

where A and C are variational parameters. This interpolative expression for u gives the 
correct asymptotic form at coincidence, r — > ("cusp conditions"). 
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u{r)\,^_^Q = const - /?^r + . . . , 



(17) 



and at very large values of r, 

«(r)U^ = -C^lnr + .... (18) 

The later expression can be easily verified for the N = 2 system. The coefficient in front of 
—In r is roughly the difference e{f3) — e(0), which is proportional to (3^ for small f3 values, 
and tends to as /? increases. 

Fig. ^ compares the Monte Carlo and Pe.sl/?) Fade results for 20 and 210 bosons as 
function of (3 in the interval < j3 < 1.3. In the computations, 10^ sweeps were used to 
compute mean values after 2 x 10^ steps for thermalisation. Note that (3 = 1.3 is well outside 
the perturbative regime. Indeed, e(1.3)/e(0) > 12 for = 20, whereas e(1.3)/e(0) > 35 for 
N = 210. 

In Fig. |I| we have also plotted (the dots) the relative difference between the Pq^^^P) Fade 
approximant and the variational Monte Carlo results. Note that in Fig. |l|(a) for = 20 
up to /3 ~ 2.5 (the solid dots) the Pq^s{P) approximant gives practically identical results 
as those obtained from the variational Monte Carlo. The relative error is less than 1.5% 
which is below the 4% level estimated previously from a comparison between successive Fade 
approximants. For larger values of /3 one approaches the Wigner limit where the particles 
form a crystal-like structure. In that limit the functional form taken for the trial function is 
expected to be no longer good and we notice that the relative error does not decrease with 
P (open dots) but stays approximately constant. Also for A^ = 210 it is shown in Fig. |I](b) 
that the relative error is less than 4%. 

ACKNOWLEDGMENTS 

AG acknowledges support from the Colombian Institute for Science and Technology 
(COLCIENCIAS). BF is an aspirant, and FMF a research director with the Flemish Science 
Foundation (FWO Vlaanderen). Fart of this work is supported by FWO, a NATO Linkage 
grant, the 'Interuniversity Foles of Attraction Frogramme - Belgian State, Frime Minister's 
Office - Federal Office for Scientific, Technical and Cultural Affairs', and the WOG on field 
theory and statistical physics. 



5 



REFERENCES 



* Electronic mail: agonzale@fisica.udea.edu. co 
° Electronic mail: bpartoen@uia.ua.ac.be 

^ Permanent address: Semiconductor Physics Institute, Gostauto 11, 2600 Vilnius, Lithua- 
nia. Electronic mail: matulis@pub.osf.lt 

* Electronic mail: peeters@uia.ua.ac.be 

^ M. H. Anderson, J. R. Ensher, M. R. Mathews, C. E. Wieman, and E. A. Cornell, Science 
269, 198 (1995); K. B. Davis, M.-O. Mewes, M. R. Andrews, N. J. van Druten, D. S. 
Durfee, D. M. Kurn, and W. Ketterle, Phys. Rev. Lett. 75, 3969 (1995); C. C. Bradley, 
C. A. Sackett, J. J. Tollett, and R. G. Hulet, Phys. Rev. Lett. 75, 1687 (1995). 

^ See, for example, S. Giorgini, L. P. Pitaevskii and S. Stringari, J. Low Temp. Phys. 109, 
309 (1997), and references therein. 

^M.-O. Mewes, M. R. Andrews, N. J. van Druten, D. M. Kurn, D. S. Durfee, C. G. 
Townsend, and W. Ketterle, Phys. Rev. Lett. 77, 988 (1996); C. J. Myatt, E. A. Burt, R. 
W. Ghrist, E. A. Cornell, and C. E. Wieman, Phys. Rev. Lett. 78, 586 (1997). 

^L. P. Pitaevskii and A. Rosch, Phys. Rev. A 55, 835 (1997). 

^Tin-Lun Ho and Michael Ma, |cond-mat/ 9703034 . 

^ V. Apaja, J. Halinen, V. Halonen, E. Krotscheck, and M. Saarela, Phys. Rev. B 55, 12925 
(1997); G. Sugiyama, C. Bowen, and B. Alder, Phys. Rev. B 46, 13042 (1992); S. Moroni, 
S. Conti, and M. P. Tosi, Phys. Rev. B 53, 9688 (1996); K. Tankeshwar, B. Tanatar, and 
M. P. Tosi, ICTP preprint IC/97/95. 
^D. R. Nelson, Phys. Rev. Lett. 60, 1973 (1989); D. R. Nelson and H. S. Seung, Phys. Rev. 
B 39, 9153 (1989); W. R. Magro and D. M. Ceperley, Phys. Rev. Lett. 73, 826 (1994); 
H. Nordborg and G. Blatter, Phys. Rev. Lett. 79, 1925 (1997). 
^A. Gonzalez, B. Partoens, and F. M. Peeters, Phys. Rev. B 56, 15740 (1997). 
9 A. Matuhs and F. M. Peeters, J. Phys.: Condens. Matter 6, 7751 (1994). 
i°A. Gonzalez, J. Phys.: Condens. Matter 9, 4643 (1997). 
" V. M. Bedanov and F. M. Peeters, Phys. Rev. B 49, 2667 (1994). 
^2 V. A. Schweigert and F. M. Peeters, Phys. Rev. B 51, 7700 (1995). 
13 A. H. MacDonald and D. R. Ritchie, Phys. Rev. B 33, 8336 (1986). 
i^S. Giorgini, L. P. Pitaevskii, and S. Stringari, Phys. Rev. Lett. 78, 3987 (1997). 
1^ D. M. Ceperley, Lectures on quantum Monte Carlo, in Spring Colleges in Computational 
Physics, ICTP, Trieste, 1997. 



6 



FIGURES 



FIG. 1. The relative differences between consecutive Fade approximants of the sequence 
PK+3,Ki(3)- We consider two systems: (a) = 20 and (b) N = 210. The dots give the rela- 
tive difference \Evmc — -P6,3(/?)|/^6,3(/5) between the best Fade approximant and the variational 
Monte Carlo (VMC) results. 

FIG. 2. Approximate scaling properties of the energy for = 90 and = 210. As a compari- 
son we also show the result for N = 20. The dots (solid dots for = 210, open dots for A^ = 90 
and stars for A^ = 20) are the results from the variational Monte Carlo calculation. 

FIG. 3. Comparison between the Pq^3{P) Fade estimates (solid curves) and the variational MC 
calculations (dashed curves) for: (a) A^ = 20 and (b) A^ = 210. 



APPENDIX: SECOND ORDER ENERGY CORRECTION FOR THE TWO 
PARTICLE INTERACTION PROBLEM 

Here we present an alternative but more accurate calculation of the coefficient C2 for 
the second order energy correction (^). Because expression (^) exhibits a rather simple 
dependence on the number of bosons we can limit ourselves to calculate the coefficient C2 
for the system of two particles with Coulomb repulsion in a parabolic confinement potential. 
The ground-state of that system will be symmetric with respect to the permutation of the 
particles what automatically takes into account the boson character of the problem. 

It is known that in this two particle interaction problem the relative and the center-of- 
mass motion can be separated. The center-of-mass {R = (ri + r2)/2) motion is the one of 
a harmonic oscillator and is easily eliminated. Consequently, we have only to consider the 
following radial Schrodinger equation for the relative motion (r = ri — r2) of particles 



with A = 13^. 

Now let us expand the eigenvalue E = Eq + XEi + ■ ■ ■ and the wave function R = 
Rq + XRi + ■ ■ ■ into powers of A. The application of the standard perturbation technique 
leads to the following expressions 



with Rq = exp(— r^/4) and the first order wave function correction obeys the following 
equation 




(Al) 



(A2) 




(A3) 
(A4) 



(A5) 
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{Ho - Eo} Ri = (1/r - E,)Ro. (A6) 
This equation has to be solved together with the boundary and orthogonahty conditions 

i?i(0) = i?i(oo) = 0, / drrRi{r)Ro{r) = 0. (A7) 

Jo 



Inserting the substitution Ri = {I]^iWri,r" + C}Ro into equation ( |A6|) we obtain the 
recurrence relation 

Wn+2 = Wnu/in + 2)^, Wi = 1, W2 = -Ei/A. (A8) 

Now iterating the above relation, withy the ortogonality condition ( |A7| ), and inserting the 
obtained result into expression (^) we obtain the second order energy correction expressed 
as a single sum of F-functions 

nr - r r(A; + l) r(A; + 3/2) 1 

-C/9 = -l / 7 



16 ^0 + l/2)r(fc + 3/2) {k + l)r(A; + 2) J 

~ 8,^„(fc + l/2)(A: + l)- ^ 

The coefficient C2 = E2/2 — ci follows from expression (|^) with N = 2. The sums in 
expression ( |A9| ) converge which enables us to obtain the coefficient C2 up to the desired 
accuracy. 
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